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T— I The energy source of the anomalous X-ray pulsars^ (AXPs) is not well understood, 
•/"J hence their designation as anomalous. Unlike binary X-ray pulsars, no companions 
^ are seen, so the energy cannot be supplied by accretion of matter from a companion 
^ star. The loss of rotational energy, which powers radio pulsars, is insufficient to power 
C^AXPs. Two models are generally considered: accretion from a large disk left over from 
^ the birth process, ^ •^ or decay of a very strong magnetic field (10^^ G) associated with a 
^ 'magnetar.'^ The lack of counterparts at other wavelengths has hampered progress, in 
• "-j our understanding of these objects. Here, we present deep optical observations of the 
^ field around 4U 0142+61, which is the brightest AXP in X-rays. The source has no 
associated supernova remnant, which, together with its spin-down time scale of ~ 10^ yr 
(ref. 5), suggests that it may be relatively old. We find an object with peculiar optical 
colours at the position of the X-ray source, and argue that it is the optical counterpart. 
The optical emission is too faint too admit the presence of a large accretion disk, but 
may be consistent with magnetospheric emission from a magnetar. 

The field around 4U 0142+61 was imaged in the V, R and I bands on 31 October 1994 using the 
Low- Resolution Imaging Spectrometer^ on the 10-m Keck I telescope on Mauna Kea, and in R and I 
on 6 September 1999 on the Keck II telescope. In 1994, the sky was not photometric and the seeing 
mediocre at 1 arcsec, but in 1999 the sky was photometric and the seeing 0.6 arcsec. Photometric 
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calibration was achieved with images obtained under photometric conditions on 23 July 2000 at the 
60-inch telescope on Palomar Mountain. 

Within the error circles of 4U 0142+61 (Figure 1), we find one faint stellar object with unusual 
colours (Figure 2). To understand its nature, and its relation to 4U 0142+61, we need to constrain 
the run of reddening with distance along the line of sight. We use two open clusters in this region of 
sky'^ - NGC 654 (at an angular offset of 17', distance of 2.7 ± 0.4 kpc and reddening Ay = 2.3 + 3.4) 
and NGC 663 (30', 2.8 ± 0.4 kpc, 2.2 + 3.1) - to infer that Ay ~ 2.7 at = 2.7 kpc. Furthermore, 
we use the bluest galaxy in our field (object G), with R — 25.6, R — I = 1.3. Because the intrinsic 
colour^ is unlikely to be bluer than {R — I)o — 0.3, the total Galactic reddening along this line of 
sight is constrained by Er^j ^ 1.1, or equivalently Ay ^ 5. 

From Figure 2, we see that most stars have brightnesses and colours consistent with those ex- 
pected for low-mass stars at distances similar to the open clusters. Star A, however, is either too 
dim or too blue. It might be bluer because it is nearer and therefore less reddened. But then it 
would be too dim to be a main sequence star and too red to be a white dwarf. The latter can be 
seen by comparison to star D, which is probably a foreground white dwarf (Figure 2). 

We proceed to estimate the temperature T and radius R of star A on the assumption that the 
object is located beyond the open clusters (i.e., d > 2.7 kpc; 2.7 ^ Ay ^ 5); below we use ^5, the 
distance normalized to 5 kpc. Bearing in mind the intrinsic variability in the shape of the reddening 
curve, ^ we find a marginal fit for T ~ 6000 K and R ~ O.llrfs Rq {Ay = 2.7), where Rq is the Sun's 
radius; and increasing ly better fits for T ~ 10^ K and R ~ 0.084 ^0 {Ay ~ 4) up to T ;^ 10^ K 
and R ~ 0.022^5 (T/IO^ K)-V2 (^Ay ~ 5.4). We note that for T ;^ lO^K the optical bands fall 
in the Rayleigh Jeans tail and the colours do not vary with T any more. Because temperature and 
reddening compensate each other, the above scaling in the Rayleigh- Jeans limit actually holds to 
within 10% for any T ^ 6000 K. 

For all temperatures, the radii are too small for star A to be a normal star within our Galaxy. 
The remaining possibihty is that star A is a reddened but hot (T ~ 5 x 10^d| K) white dwarf. This 
is extremely unhkely unless it is actually 4U 0142+61 itself, as we discuss below. 

From the above, we conclude that in all likelihood star A is the optical counterpart to 4U 0142+61. 
Below, we discuss our observations in the framework of the models that have been proposed for AXPs. 
The optical extinction Ay to 4U 0142+61 is between 2.7 and 5.4 mag, as estimated from the X-ray 
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dust scattering halo and absorption column density, respectively^^'^^ (where the latter is hkely the 
more reliable^"); therefore, the source must be at a distance exceeding ~2.7kpc. 

We first consider the possibility that AXPs are isolated neutron stars accreting matter from a 
disk, presumably composed of supernova debris. ^'^ In this model, the optical emission arises in the 
accretion disk, mostly due to reprocessing of the X-ray irradiation^^ (Figure 3). The isotropic X-ray 
luminosity of 4U 0142+61 is Lx = Awd^ f^""^^^ = 1.8 x lO^^dlergs'^, where /^"^^^^ is the flux in the 
energy range 0.5-10 keV after correction for interstellar absorption (see Table 1). Wc can estimate 
the run of temperature with disk radius^^ r by T(r) ~ 5OOO(//O.25)^/^(i5^'^(r/i?0)~^/'' K. Here, the 
factor / parametrizes uncertainties in vertical disk structure and the fraction of impinging X-ray 
emission reflected by the disk. 

The flux at a given frequency and for given inclination is obtained by integrating the emission from 
the accretion disk (assumed to be optically thick) from an inner radius, Tin, to an outer radius, rout; 
see ref. 12 and Figure 3. For a simple estimate, we use that the optical flux will arise predominantly at 
those radii at which the emission peaks in the optical band, i.e., where T ~ 5000 K. For 4U 0142+61, 
this would be at r ~ l{f /0.25)'^^'^dt/^ Rq, which is much larger than the limit of i? ~ 0.11^5 i?© 
derived above, and leads to a brightness much larger than observed (see Figure 3). Thus, this model 
is excluded. 

The optical emission from the disk can be reduced by appealing to a larger Tin or a smaller rout- 
The former, suggested^^ in the context of earlier limits, would result in extremely red emission 
(because of the absence of the hotter inner region). This is incompatible with the observations. 

A reduced rout would follow naturally in a model less often considered, in which the AXP are 
compact binaries. We find that our observations require rout — 0-05dl°^^^ (f /0.25)~'^^^^ Rq] see 
Figure 3. This is very small, but not unprecedented: the X-ray binary with the tightest known 
orbit, 4U 1820—30 (an 11 min orbital period) has a similarly small radius and similarly high ratio 
of X-ray to optical flux; see Table 1. However, the X-ray spectrum of 4U 0142+61 is rather different 
from that of 4U 1820—30 and other accreting sources. Furthermore, such compact binaries are not 
expected to be associated preferentially with supernova remnants. Thus, we consider this model 
unlikely. 

As mentioned, the optical data are consistent with a hot white dwarf. This would be expected in 
another model rarely considered, in which AXPs are massive and magnetized (10^ G) white dwarfs 
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rotating on the verge of breakup, presumably formed in the merger of two ordinary, approximately 
0.6 white dwarfs. In this model, the ultimate source of energy, as in radio pulsars, is rotation. 
For 4U 0142+61, the inferred rotational energy loss, — -Erot = 47r^/P/P^ ^ lO^^ergs"^, is clearly 
sufficient to account for the X-ray luminosity; here, / — kMB? ~ 10^°gcm^ is the moment of inertia 
appropriate for a hot white dwarf with mass M ~ 1.3 Mq, radius R ~ 0.007 Rq and gyration constant 
k ~ 0.14. P is the spin period and P its derivative. If some fraction /3 of E^ot goes into heating, as 
might result if the white dwarf rotated differentially, we expect a surface luminosity L = AnR'^aT^ — 
-pErot, and thus a surface temperature T ~ 4 x 105(/?/0.5)^/^(A;/0.14)^/^(M/1.3 Mo)^^ k (we note 
that T does not depend on it!). This is sufficient to produce the optical emission for a source at 

— 0.6{R/ 0.007 Rq) (see Figure 3). We note that a possible descendant may already have been 
identified^^: the 1.3 Mq white dwarf RE J0317-853 with P = 725 s, 5 ~ 5 x 10^ G, and a cooling age 
of a few 10* yr. In this model, however, what mechanism may cause the X-ray spectrum is unclear, 
the thermal emission being far too soft (see Figure 3). Furthermore, the association of other AXPs 
with supernova remnants again seems puzzling. 

This leaves us with the magnctar modcl,'^ in which both the power law component of the X-ray 
emission and the optical emission would be of magnetospheric origin. Unfortunately, there are no 
detailed models. Radio pulsars with X-ray and optical detections have rather smaller fx/ f opt (see 
Table 1), but their X-ray spectra are unlike those of AXPs. From Figure 3, it appears that the 
simplest spectral energy distribution would be obtained if the extinction were near the high end 
{Ay ^ 5), so that the emission could peak at a few 10^^ Hz (~ 100 eV) and be like a power law 
fi, oc with index a ~ 2 (or even 2.5) in the optical. This could arise if the magnetospheric 
emission is self-absorbed at optical frequencies. A possible problem for the magnetar model is the 
lack of X-ray emission from the one radio pulsar with similar P and P (ref. 19). However, it may 
well be that, P and P are not reliable measures of the strength of magnetic fields for AXPs and that 
other properties determine whether a source is an AXP. 

In spite of the current uncertainties, the optical identification offer us new insights into these 
enigmatic objects and motivates new observations and searches for optical pulsations. 
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Table 1. Comparison with other high energy objects 

Type of object: AXP: anomalous X-ray pulsar; PSR: radio pulsar; XRB: X-ray binary; INS: isolated 
neutron star. 

Although RX J0720.4— 3125 has a similar pulse period and X-ray to optical flux ratio, its X-ray 
spectrum is very different from those of the AXPs. 

The observational data have been collected from the literature. A full set of references can be found 
on http:/ /www.astro.uu.nl/~hulleman/u0142. The unabsorbed X-ray fluxes are for the energy range 
of 0.5-10 keV. They have been inferred from the best-fit spectral model with the help of the PIMMS 
software package. For 4U 1626—67, constant pulsed flux is assumed in the estimate of {R — /)o, and 
for 4U 1820—30, we assumed V — R — to estimate R. For the conversion to fluxes, we used the 
BesseU zero points,^^ where V ^0 corresponds to 3600 Jy, = to 3060 Jy, and 7 = to 2420 Jy. 
For reddening corrections, we used {Ay^^, Ar, Aj) = (1.015, 0.819, 0.594)^4^, where Ay^^, Ar, and Aj 
are reddenings in the Landolt V, R, and 1 fllters, and Ay is the reddening in the Johnson V band; 
see Appendix C in ref. 21 for a discussion. 

For reference, the flnal column lists estimated distances; the colons indicate the uncertainty, zero 
colons indicating a factor up to 1.5 uncertainty, one colon a factor two, two colons a factor three. 
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Figure 1. Keck images around the X-ray position of 4U 0142+61. The large panel shows an overview 
of the R-band image. Our candidate, A, as well as two other blue objects, D and G, are indicated 
with labels just to the right of their images. The small bottom panels are close-ups of the 1-band 
images of A and G, which show that G is extended and probably a galaxy, while A is not. The panels 
on the right show the V, R, and I images around star A. The error circle indicates the X-ray position 
derived from analysis of Einstein HRI observations (aj2ooo = 0lM6"'22!03, 5j2ooo = +61°45'04('3, 95% 
confidence level radius of 3.9 arcsec; ref. 36). The Einstein positions have proven to be highly reliable, 
but for comparison we display in the I-band panel also the error circles derived from observations 
with the EXOSAT LETT (end figures 22?40 and 10{'2; 4.8 arcsec radius; ref. 36) and ROSAT PSPC 
(21?93, 44'57^'7; 6 arcsec radius). The ROSAT position of 4U 0142+61 was derived from archival data, 
using a boresight correction inferred from two other X-ray sources in the field, for which the optical 
counterparts are known^^'^^ (LS I +61° 235 and LP 80-77; our technique is described in ref. 14). 
The astrometry was done as in ref. 14, by measuring 107 stars from the USNO-A2.0 catalogue^^ on 
a short R-band exposure (root-mean-square residuals of 0.26 arcsec in each coordinate, consistent 
with the USNO-A2.0 measurement errors). From the 1-band image, we infer a position of star A of 
a;j2ooo = 01^46"'22!41, (5j2ooo = +61°45'03'.'2. We estimate that this position is on the International 
Celestial Reference Frame to about 0.2 arcsec. The other images give consistent positions. Our 2a 
upper limit to the proper motion is 0.03 arcsec yr~^. 
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Figure 2. Colour- magnitude and colour-colour diagrams for stars near the position of 4U 0142-1-61. 
For comparison, the magnitudes and colours expected for stars^*^ at the distance, reddening, and age 
of the nearby open cluster NGC 654 (2.7 ± 0.4 kpc, 2.7 mag, 50Myr; rcf. 7; solid line) are shown, 
as well as a cooling track for a 0.6 white dwarf'^^ at the same reddening and distance (dashed 
line). The arrows indicate the effect of increasing the reddening by A Ay — 1. All optical images 
were corrected for sensitivity variations using dome flats, and a small correction for non-linearity 
of the detector was applied to the 1999 observations. For the 1999 I-band observations, a 'fringe' 
frame was constructed by taking the second-lowest value for each pixel from seven dithered images, 
which was subtracted from the individual images to remove interference patterns. Instrumental 
magnitudes were measured using the DAOPHOT II package^^ on the 1994 V band and the 1999 R 
and I-band images, using only those regions not affected by scattered light from bright, overexposed 
stars. The 1994 V-band images were taken at a different pointing and rotation angle, which caused 
different parts on the sky to be affected by bleed trails, and so on. Therefore, we obtained V-band 
measurements only for stars in the central part displayed in Figure 1. The photometric calibration 
was done relative to the 60-inch data, which were calibrated in turn using exposures of the standard 
fields PG 1657—042 and NGC 7790.^^ The R and 1 calibration was verified using Keck images of the 
standard field PG 0231-1-051.^^ We estimate the uncertainties in the zero points of the magnitudes 
and colours to be about 0.03 mag. For star A, we find R = 24.99 ± 0.07, V - R ^ 0.63 ± 0.11, 
i? — / = 1.15 ± 0.09. By comparison with the 1994 images, we find that its brightness was constant 
to within 0.2 mag {2-a) in R. 
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Figure 3. Energy distribution for 4U 0142+61. At low frequencies (lO^'^-lO^^ Hz), the points marked 
V, R, and I indicate the observed V, R, and I-band fluxes. The vertical error bars reflect the uncer- 
tainties, while the horizontal ones indicate the filter bandwidths. The set of points above the mea- 
surements indicate dcreddened fluxes for Ay = 5.4, as inferred from the X-ray column density. 
The errors include a 3% uncertainty in the reddening correction.^ At high frequencies (10^^-10^® Hz), 
the crosses show the incident X-ray spectrum as inferred from ASCA measurements.^'^ The diamonds 
show the spectrum after correction for interstellar absorption, and the two thick dashed curves show 
the two components used in the fit^*^: a power law of the form F,, = 103(/;,z//l keV)"^'^'' yuJy and a 
black body with T = 4.4 x 10^ K and R = 12d^ km. The latter component is extrapolated to lower 
frequencies to show it cannot reproduce the optical fluxes. Drawn in thin solid lines are models for 
the optical emission. The curve marked disk is for an accretion disk that has inner radius equal 
to the corotation radius i?co = {GMP"^ /Att'^Y/^ = 0.010 i?© (for M = 1.4 M©, P = 8.7s), outer 
radius of 10^^ cm, and inclination of 60 degrees. Both irradiation and viscous heating are taken into 
account the former dominates in the optical and the latter causes the bump in the ultra-violet. 
The calculation is for d = 5 kpc, but since the disk luminosity scales almost linearly with the X-ray 
luminosity, the result is not sensitive to distance. The optical fluxes expected in the accretion model 
are greatly in excess of those observed. Only for a truncated disk with a small outer disk radius, 
'^out ^ 0.1o?5°^^^(//0.25)~^/^^ i?0, can the observed fluxes be reproduced, as is shown by the dashed 
curve. The thin solid curves marked BB, H, and He are black body, pure hydrogen and pure helium 
white dwarf model atmospheres for T = 4 x lO^K. The normalization to the optical data implies 
R = 0.011^5, O.Ond^ and 0.015^5 i?0, respectively. None of the spectra can reproduce the X-ray 
emission. 
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